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The Brownian dynamics of a single microorganism coupled by chemotaxis to a diffusing concen- 
tration field which is secreted by the microorganism itself is studied by computer simulations in 
spatial dimensions d = 1,2, 3. Both cases of a chemoattractant and a chemorepellent are discussed. 
For a chemoattractant, we find a transient dynamical arrest until the microorganism diffuses for long 
times. For a chemorepellent, there is a transient ballistic motion in all dimensions and a long-time 
diffusion. These results are interpreted with the help of a theoretical analysis. 
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I. INTRODUCTION 

Chemotaxis [H HI [S] and Brownian motion [H |51 |S] be- 
long to the key processes which govern the motility of 
microorganisms (e.g. bacteria, amoeba and endothelial 
cells) [7]. In the simplest approach, the microbe "smells" 
a chemical and moves along the gradient of the concen- 
tration field of the chemoattractant in order to reach effi- 
ciently the secretion source of the chemical. The opposite 
case of negative chemotaxis is realized in case the mi- 
crobe intends to avoid another object which is secreting 
the chemical [8]. This chemotactic drift is superimposed 
to stochastic motion due to fluctuations associated with 
the active process of self-propulsion of the microorganism 
P'. Chemotaxis can lead to clusters of aggregated bac- 
teria [10, 111 which are still emitting chemoattractant. 

Here we study the self-coupled situation where the mi- 
croorganism "smells" itself, i.e. it reacts chemotactically 
to its own secreted chemical. This "autochemotaxis" is 
realized in aggregated clusters of different bacteria if the 
aggregate is considered as a net particle. Another real- 
ization is a single bacterium which has both an emitter 
and a sensor of the same chemical. Tsori and de Gennes 
[12j have studied a simple model for this situation in dif- 
ferent spatial dimensions d and find self-trapping of the 
bacterium in its own chemoattractant cloud for d — 1,2 
but not for d = 3. This means that for low dimension- 
ality the bacterium is fooled by its own secretion such 
that it is getting localized for long times. In a subse- 
quent numerical study of a model microbe coupled to its 
own chemoattractant secreted at constant rate, Grima 
[131 [14] calculated the long-time dynamics in various di- 
mensions d and found long-time diffusive behavior even 
for c? = 1 , 2 in contradiction to Ref. [T2] . Grima also stud- 
ied the case of negative chemotaxis and finds in all dimen- 
sions long-time diffusion or ballistic motion depending on 
the strength A < which couples the microbe's driving 
force to the gradient of the chemical concentration field. 
Grima predicts that if |A| exceeds a critical value Ac, then 
the long-time motion is super-diffusive. 

In this paper we revisit autochemotaxis by studying a 
model which is similar but not identical to that proposed 
by Grima [I3j. In the model of Grima, a finite global 



extinction rate (A) of the secreted chemical is present. 
Here we focus on the case A = 0, which can be realized 
in experiments (see e.g. 0). By extensive computer sim- 
ulations, we study different spatial dimensions d = 1, 2, 3 
and both cases of positive and negative chemotaxis. Our 
results are summarized as follows: Consistent with Grima 
|13j . we find for positive chemotaxis (i.e. for a coupling 
parameter A > 0) long-time diffusive motion. In addition 
we find dynamical transients before reaching the long- 
time diffusive limit. During the transients, the dynam- 
ics of the microorganism is strongly reduced resulting in 
almost dynamical arrest. The crossover time from inter- 
mediate arrest to long-time diffusion grows strongly with 
the coupling A. Therefore the idealized analysis of Tsori 
and de Gennes who predicted self-trapping (i.e. a com- 
plete dynamical arrest) is manifested by long transients 
for very strong couplings |15] . The averaged mean-square 
displacement as a function of time does not exhibit a uni- 
versal slope in this transient regime, but the actual mean 
slope is decreasing with an increase in the coupling A. 
The transient behavior is most pronounced in one di- 
mension but weakened considerably in three dimensions. 

For negative chemotaxis (i.e. dynamical self-avoiding 
of the microbe), on the other hand, we find transient bal- 
listic motion in all dimensions. In d = 2,3 we observe a 
long-time diffusion for all coupling strength. Such long- 
time diffusive motion in d = 1 , though not directly ob- 
served in the simulation, is possible by finite probability 
of changing the direction of motion (left to right) at long 
time scales for non-zero temperature. According to Refs. 
[13], the critical coupling Ac above which ballistic long- 
time behavior is found depends on the global extinction 
rate A but stays finite when A ^ 0. One reason for the 
discrepancy is because noise has not been completely in- 
cluded in the earlier treatments |13j . while solving for 
the integrals concerning the non-Markovian chemotactic 
force. We take note of the effect of noise on the mi- 
crobe's trajectory in an appropriate place in this paper 
(for another example demonstrating the importance of 
noise, see Ref. [9j). Again we address the transient be- 
havior and find an intermediate time window where the 
super-diffusive motion is found between a short-time and 
long-time diffusive behavior. This motion we observe is 
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similar to persistent random walk of a microbe which can 
be mapped to the worm-like chain model [71 [16] . In the 
particular case we consider, the persistence length is seen 
to depend on the microbe's coupling with the repellent. 

Our predictions can in principle be verified in experi- 
ments on aggregates of bacteria. For many bacteria our 
model reduces to particles interacting via gravitation- 
kind potentials for d = 3. Therefore our analysis might 
have applications for Brownian dynamics of gravitational 
matter [171 [THl . Further generalizations of our model 
are to predator and prey models possibly leading to in- 
teresting spatiotemporal delay effects, see e.g. pOj . 

Our paper is organized as follows: in section II, we 
propose the model of a microorganism coupled to its own 
chemoattractant/chemorepellent, provide the simulation 
details and point out experimental situations to compare 
typical estimates of parameters used in the model. In 
section III we present the results of our investigation. 
In section IV we explain our findings with simple the- 
oretical analysis. We conclude the paper in section V 
by discussing the main points of our findings, compar- 
isons of our results with relevant practical cases of self- 
propulsion, and future directions of our research. 



II. THE MODEL AND SIMULATION DETAILS 

Model: The microorganism is modeled as a point par- 
ticle which undergoes completely overdamped Brownian 
motion with an effective temperature in a medium 
with viscosity coefficient 7. This 'particle' is assumed 
to emit a chemical, with which it is self-coupled, con- 
tinuously in time. In this context, the word 'particle' is 
thereby taken to represent the chemotactic agent of in- 
terest - an idealized microorganism coupled with a self- 
emitted chemical field. Fluctuations associated with the 
active process of self-propulsion of the microbe is mod- 
eled by the effective temperature parameter /3~^ in our 
system. The time evolution of the density field p{r,t) 
of the continuously emitted chemical is thus governed by 
a diffusion equation with a source term which depends 
upon the instantaneous position ri,(t) of the moving mi- 
crobe: 



dp{r,t) 
dt 



(1) 



Here, the constants Xe and Dc are the rate of emission of 
the chemical and the diffusion constant of the chemical 
in the medium, respectively. 

In the absence of chemical, the microbe diffuses non- 
chemotactically in the medium with an effective free dif- 
fusion constant D = 1/(7/?). However, with the presence 
of the emitted chemical the resulting 'chemotactic' be- 
havior depends on the nature of the self-coupling of the 
microorganism to its chemical field, i.e. whether it moves 
'up' or 'down' the chemical density gradient. We study 
both cases by simply modeling the self-coupling 'force' to 
be proportional to the gradient of the chemical density 



field Vp{r,t), and the proportionality constant A deter- 
mines the strength as well as the nature of the coupling. 
In reality, of course, chemotaxis can be more complex 
involving temporally sampling of the concentration field 
and a biased random walk [71[1TJ[52]. In this context, we 
note that the 'chemotactic force' imitates the net effect 
of chemotactic movements on a phenomenological level. 
The details of the actual propulsion mechanism [25] of 
the chemotactic agent and the effect of the solvent flow 
field on the diffusing chemical are not taken into account. 
Within our simple model, positive and negative A natu- 
rally generate the cases of positive and negative chemo- 
taxes, respectively. We are thus lead to the following 
idealized model of a chemotactic agent 



7r6(i) =F(rfc,^)+T7(<). 



(2) 



Here, ri{t) is an effective noise specified by {ri{t)) = 
and {'ni{t)rij{t')) = 2jf}~^5ijS{t — t'), with i and j refer- 
ring to the spatial components of the noise vector. This 
noise term is assumed to effectively take care of all non- 
equilibrium fluctuations that may be associated with the 
active process [2T, of self-propulsion, in absence of the 
chemical. F(rb,t) denotes the model 'chemotactic force' 
taken to imitate the systematics of the effective chemo- 
tactic movement of the microbe at the position r?, at time 
t due to the chemical secreted all along the trajectory 
traversed in the past. It is obtained by analytically solv- 
ing Eq. ([T| for the density field p{r,t) by the method of 
Green's function, and subsequently calculating the gra- 
dient Wp{r,t). The 'force' at a time instant is dependent 
on the entire previous path history of the microorganism, 
thereby generating a strongly non-Markovian dynamics. 
However, owing to a physical memory time (to) associ- 
ated with the microbe to sense its chemical, the part of 
the trajectory in this most recent time tQ, i.e. for all ri,{t') 
with t — to < t' < t, does not contribute. The physical 
import of this is that there is a finite time delay to , how- 
ever small, between the act of secreting chemical by the 
microorganism and the act of responding to it, during 
which the sensor gets to activate. With the introduction 
of the memory time, to, the 'chemotactic force' at time t 
and at position r becomes 



F{r,t) = -2AAe 



4D^\t-f\ {iTTDc\t-V\y/^'^ 



(3) 

where d is the dimensionality of the embedding space. 
Evidently form Eq. ^ , for the mathematical case of to — 
0, the 'force' becomes divergent. 

Simulation details: We performed extensive Brownian 
dynamics simulation [26] for this non-Markovian process 
of a microorganism moving by autochemotaxis. We mea- 
sured time in units of tq = A~^, all lengths in units of 
lo — (-s/UD^/Ae)^/^ and energies in units of /3~^. The 
coupling strength A is measured in units of /?~^Zo'^- Thus, 
we set Ae = 1, Zq = 1 f^^d (3 = 1 for convenience. 
Further we considered the physical situation when the 
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microorganism diffuses at a much slower rate compared 
to the emitted chemical in the medium [271 EH] , taking 
D — 0.1 Iq' Itq and fixing the ratio D/Dc = 0.01. In our 
Brownian dynamics simulations we used to — 0.001 tq. 
The Langevin equation [Eq. Q] is solved with a discrete 
time step A< = 0.0001 tq. Space is, however, continuous. 

Connection to experiments: In order to get an estimate 
for the coupling strength A in our units, we note that 
the typical value of the ejection rate of chemical from 
a microorganism is Ae - 10^ molecules/s [El [29], and 
usually D/Dc ^ 10"^ - 10"^ HZl [55]. In three spatial 
dimensions, for example, the chemotaxis of Dictyostelium 
to shallow cAMP gradients [301 EI] with typical values of 
Vp - 0.01 nM/fim, ~ 300 ^im^/s, D/Dc - 10"^ 
and moving with steady state velocity v 0.2 fim/s, 
yields A - 10-* 

Chemotactic Microglia cells [28 move at a speed of 
V ^ 2/im/min in a spatial gradient Vp ~ 0.003 nM / 
of a chemoattractant (IL — 1/3) which is secreted at a 
rate of Ae ~ 200 molecules/min and diffuses with Dc = 
900 /im^ /min. The coupling strength in this case is A ~ 
10 l3~^lo for an effective non-chemotactic diffusion con- 
stant Z? ~ 33 pm^/min of microglial cells, due to random 
motility in the tissue. The chemoattractant has a low de- 
cay rate of A 0.003 — 0.03 min^^. Microglial cells are 
also known to respond to a chemorepellent (TNF — a) 
they produce, with similar production, diffusion and de- 
cay rates as the chemoattractant. Our estimate of the 
chemoattractant gradient was based on the value of the 
chemotactic coefficient ^ 780 /im^ nM~^ min^^, a ra- 
tio between cell velocity and chemical gradient, used in 
Ref. (28J. The corresponding value for the repellent is 
not known. 

Again, for E. coli of size ^ 1 /xm, swimming at 
f ~ 20 fj,m/s in the spatial gradient Vp ^ 0.1 fj,M/^m 
of a chemoattractant diffusing with Dc — 10^^ cm'^/s, 
in a medium of viscosity 10"'^ Pa.s [32l|33]; the coupling 
strength is A ~ 10"! f3-^lo^. The non-chemotactic diffu- 
sion coefficient of the bacterium is D = 6.6 x 10"^ crv? j s. 
The time required by the chemical in this case, to diffuse 
a length equal to the size of the bacterium, is of the order 
of 0.1 Tq. The memory time, needed by the bacterium 
to respond to the chemical stimulus can be much smaller 
than this time. 

In all the above calculations, the effective non- 
chemotactic diffusivity D was used to express the energy 
unit fi~^ = 7-D. 



III. RESULTS 

We now investigate the nature of the dynamics in all 
the dimensions and for both cases of positive and negative 
A, examining the model microorganism from some initial 
reference point taken as the origin, i.e. V},(t = 0) = 0. For 
this purpose we computed the mean square displacement 
of the microbe as a function of time and averaged over 10'^ 
realizations for each case. We checked that the system is 



in a steady state, and have also performed a steady state 
averaging of the mean square displacement. We illustrate 
our findings below. 



A. Positive autochemotaxis 

Upon examining the motion of the microbe in a 
chemoattractant (A > 0) in one dimension (d = 1), 
we found signatures of long-time diffusion (((rf,(t) — 
rh(0))^) ^ t) with a modified diffusion constant Di. The 
value of Di depends on the strength of the coupling A, 
and decreases with increasing A. For very high A values it 
is computationally difficult to obtain the long time diffu- 
sive behavior; but we obtained an upper estimate of the 
diffusivity for lower A values from a fit to the obtained 
data. In the opposite limit, i.e. at very short times, 
the microbe's dynamics is also diffusive with the non- 
chemotactic diffusion constant D. The departure from 
this behavior occurs at times dependent on the coupling 
strength: the stronger the coupling, the dynamics be- 
come markedly history dependent and the microorganism 
deviates from the non-chemotactic diffusion faster. At 
intermediate times we find long period of sub-diffusive 
crossover regime, showing signatures of a transient dy- 
namical arrest at high A-values. The crossover time also 
increases with the coupling strength. Fig. [ija) shows 
the mean square displacement as a function of time in 
d = 1 , for three values of A and compared with the non- 
chemotactic diffusion. 

In (i = 2, we find similar long-time diffusive behavior. 
The crossover times from early-time non-chemotactic dif- 
fusion to the long-time modified diffusion, is greatly re- 
duced for a given A-value as compared to the one di- 
mensional case. This feature is attributed to the effect 
of fluctuations. The crossover time to the final long-time 
diffusion, however, increases with coupling strength. The 
intermediate sub-diffusive regime is also diminished in 
time. The mean square displacement as a function of 
time, in this case, is shown in Fig.jljb). 

In Fig. [ijc) we show the case for d = 3, and in consis- 
tence with our expectation we find similar long-time dif- 
fusion with a further reduced crossover time as compared 
to the low-dimensional cases. Thus the transient dynam- 
ical arrest becomes more prominent at lower dimensions 
and at stronger couplings. The A-dependence of the mod- 
ified diffusion constant Z?;, is shown in Fig.[ljd), for all di- 
mensions, relative to the corresponding non-chemotactic 
diffusion constant D. 



B. Negative autochemotaxis 

For the repulsive case, when the model microorgan- 
ism gets repelled by its ejected chemical (A < 0), we 
found a short-time non-chemotactic diffusive motion with 
a crossover to a ballistic behavior (((r{,(t) — rb(0))^) ~ t^) 
in one dimension, for all values of the coupling strength 
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FIG. 1: (color online). Mean square displacement {{rb{t) — rb(0))^) of the microorganism as a function of time t with chemoat- 
tractant in (a) d = 1 with A = 20, 500, 5000; (b) d = 2 with A = 1000, 10000, 40000; (c) d = 3 with A = 1000, 10000, 40000. The 
non-chemotactic diffusion reference lines are also indicated as 2Dt, ADt and 6Dt correspondingly for d = 1, 2, 3. Reference lines 
(thick dotted) are used to indicate the long-time diffusive behavior (~ t) wherever possible. The relative long-time diffusivity 
Di/D is shown as a function of A in (d) for d = 1,2, 3. The reference line (thick dotted) shows a power law scaling behavior 
~ 1/A^ (see text). 



(Fig. [2|a)). The time of commencement of the ballistic 
behavior, however, depends on A. For weak coupling (i.e. 
low |A|), the dynamics resembles the non-chemotactic dif- 
fusion for longer times before finally going over to the 
ballistic dynamics. The velocity of this ballistic motion 
is given by the time derivative of the root mean square 
displacement, ^ yj ii^bit) — rh{0))^), and this velocity is 
seen to increase with increase of the coupling strength 
|A|. We argue that the microorganism can change the 
direction of its motion in one dimension due to a non- 
vanishing finite barrier crossing probability. Under such 
circumstance the motion will be diffusive at very long 
times (not seen in the simulation), with higher diffusion 
constant. 

In higher dimensions, d — 2 (Fig. [2j^b)) and d — 
3 (Fig. [2]jc)), we observed a crossover from the non- 
chemotactic diffusion to a long-time diffusion for all cou- 
pling strengths, with a ballistic transient dynamics at 
intermediate times for high | A | . The velocity of the tran- 
sient ballistic motion increase with increasing coupling 
strength. The time duration of this ballistic transient 
as well as the modified long-time diffusion constant are 



also A-dependent, both increase with increase in |A|. In 
Fig.|2jd) we show this modified diffusion constant for the 
repulsive case relative to the non-chemotactic diffusion 
constant as a function of |A|, in d = 2, 3. 



IV. THEORY 

Our findings can be understood qualitatively and 
partly quantitatively by relatively simple theoretical con- 
siderations, which are presented in the current section. In 
subsection |IV A( a scaling law for the diffusion constant 
for the case of positive chemotaxis, Di cx D^^^Iq^^ / , 
is derived based on a simple rate theory. The same scal- 
ing Di (X had been predicted for the slightly differ- 
ent model , which includes evaporation, by Newman and 
Grima [3¥ and Grima [13] • In subsection IV B we present 
a theory, which quantitatively predicts the long-time dif- 
fusion constant in the case of negative chemotaxis for all 
coupling strengths. 

Both approaches, for negative and positive chemotaxis. 
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FIG. 2: (color online). Mean square displacement {{rb(t) — rb(0))^) of the microorganism as a function of time t with 
chemorepellent in (a) d = 1 with A = -100, -300, -50000; (b) d = 2 with A = -1000, -40000, -100000; (c) d = 3 with A = 
— 1000, —40000, —100000. The non-chemotactic diffusion reference lines are also indicated as 2Dt, ADt and 6Dt correspondingly 
for d = 1,2,3. Reference lines (thick dotted) indicating the ballistic (~ t^) and the long-time diffusive (~ t) dynamics are 
shown as guide to the eye. The relative long-time diffusivity Di/D is shown as a function of |A| in (d) for d = 2,3. The points 
represents the actual data obtained from simulations, the lines correspond to a semiquantitative theory (see text). 



attribute the long-time diffusion to small fluctuations 
about the respective steady states in the case of zero 
noise. In each subsection we therefore first present the so- 
lutions to the equation of motion (|2| without fluctuations 
(i.e., at zero temperature or inflnite coupling strength A), 
before discussing the influence of small fluctuations on 
the long-time dynamics. 



A. Positive autochemotaxis 

In the case of strong positive autochemotaxis the 
model microorganism is trapped within its own secretion, 
which effectively provides an attractive external potential 
at the microbe's location rb(t). For zero noise, the mi- 
crobe is at rest, i.e., rh{t) — Vh, and does not experience 
any force. If, on the contrary, the microbe was at time t 
instantaneously placed a distance r — Tf, away from the 
location, which it occupied at all earlier times t' < t, it 
would feel a force Fs(r), which is obtained analytically 
by evaluating Eq. ^ (see also Fig. |3|: 




X 



FIG. 3: (color online). The force |Fi,(r)|, which a microor- 
ganism would experience if it was at time t instantaneously 
placed a distance r — r;, away from the location, which it oc- 
cupied at all earlier times t' < t, plotted as a function of 
X — \r — rb\/{2\/Dcto) for all dimensions d — 1,2,3. Rep- 
resentative parameters chosen are Dc = 0.1, to = 0.01, and 
A = 10^ 
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'2D. 



F.(x) = 



erf {x) 

ATTDlta \ 
X 



^erf(x) 



d=l 

d=2 

d=3. 



(4) 

Here, x — (r — r^) / {2^/Dcto) is the dimensionless distance 
from the former position r;, of the microbe and x = |x| is 
its absolute value. Locally, i.e., for distances |r — r^l <C 
2y/DctQ, the force is linear in the distance and given by 



|F.(r)| 



A|(r-rfc)|, 
{87TD%y^ X\{r- 



d=l 
d=2 
d= 3. 



(5) 



Without fluctuations, the microbe never experiences 
the described force field. However, if the noise is nonzero 
but small, the microbe eventually walks up the locally 
parabolic walls of its self-generated potential for a short 
time T, which, in turn, might lead to a shift of the po- 
sition of the minimum Tq of the self-generated attractive 
potential by a distance Arg = ro(i -I- r) — ro(i). Clearly, 
the minimum is only displaced, if the duration of the ex- 
cursion T is larger than the memory time to and if the 
microbe excurses predominantly in one direction. After 
time t + T the microbe position might relax to the new 
minimum, rf,(t' > t + t) ^ ro{t + t). By such a fluctu- 
ation, the microorganism effectively manages to move a 
distance Atq within the time r. 

Most relevant for the long-time diffusion are those fluc- 
tuations, which yield a large displacement Aro and still 
occur at a high rate 77? < 1/t. These fluctuations are 
regarded to constitute the relevant mean steps in an ef- 
fective continuous time random walk with the desired 
diffusion constant Di ^ Atq/t. 

In this subsection, we only attempt to obtain a scal- 
ing law for the long-time diffusivity. We therefore re- 
strict our consideration to a subset of very simple fluc- 
tuations, which are believed to be representative for all 
fluctuations relevant for diffusion. In particular, we con- 
sider a "jump" process at time t — 0: At all earlier 
times, t < 0, the microbe is resting at the location of 
the self-generated potential's minimum at the origin, i.e., 
rh{t < 0) = ro(0) = 0, before it undergoes an excursion 
to a new location rb(e < t < t) — Arf, within a time e; 
the latter timescale is assumed to be small with respect 
to the residence time r. During the latter time span r 
the microbe stays at the new position, where it resists the 
force due to the secretion from earlier times. In the case 
of a large coupling strength A the process of getting to 
the new location during the time span of duration e is ir- 
relevant. After time r the position of the model microbe 
is deterministically relaxing towards the new minimum 
of the chemoattractant at Aro(T, Ar^) < Ar?,, which is a 
function of the parameters r and Ar^, only. Summariz- 
ing, the simple pathway is described by 



rb(0)=O, t<0 
Tbit) = <( Arfc , e<t<T 
Aro(r), t>T. 



(6) 



The probability of the described fluctuation is propor- 
tional to the Arrhenius factor RS] 



p [t, An] cx exp [-S{t, Arfc)/4] , 



(7) 



with Arfc — |Arf,| and ^(r, Arf,) the Onsager-Machlup 
action 



5(r,Ar6) =7 / dt\h{t) - j~^F{t)\ 
Jo 



(8) 



Neglecting the initial process of moving the distance Arj, 
and also ignoring the relaxation of the chemoattractant 
during the time r the action is approximately given by 



S{t, An)^j-^T[Fs{Arb 



which, together with Eqs. Q and ([5| yields 
S{T,Arb) oc 



An < 2^/DJq 
An > 2y^D^ . 



(9) 



(10) 



According to Eq. ([t]), a minimum requirement for the 
described fluctuation to occur frequently is that S{T,An) 
does not exceed a value of the order of 1, i.e., S{t, Ar^) < 
1. As Eq. ( 10 1 is a strictly monotonically increasing func- 



tion of T and An, this constraint is equivalent to the 
equality 



S{T,An) = 1. 



(11) 



Diffusion is believed to be governed by those random 
displacements of the microorganism, which maximize 
the shift of the potential's minimum Aro(T, Ar^) = 
|Aro(T, Ar;,)| subject to this constraint. 
We will shortly see that for large coupling constants 
> i:'^+4o"^7 small excursions Arf, < 2^/25^ are 
most relevant. In this latter limit, Aro is obtained ana- 
lytically as 



Aro(T,Arf,)~ 1 - (io/r) 



^d/2 



Arf,, 



(12) 



The maximum of Arg subject to the constraint of 
Eq. ( 11 ) is given by 



Ar 



1 



-A7 



b J 



(13) 



at an optimum displacement of the microbe and an op- 
timum residence time 



Art. 



btn 



(14) 
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with prefactors a = 2-\ 3'^^^, 2-^/3 and 6 = 4, 3, 2^^^, 
in one, two, and three dimensions, respectively. Clearly, 



Eq. ( 14) fulfills the before mentioned assumption of small 
excursions m the limit of large > D'^+H^~^j. Eq. ^ 



yields the desired scaling behavior for the long-time dif- 
fusivity Di cc Ar^/r, 



Di cx 



(15) 



In conclusion, the simple theory predicts an inverse 
quadratic dependence of Di on A in the strong-coupling 
limit, which has also been observed in the computer simu- 
lations [see Fig.[T|d)]. The same scaling had already been 
predicted for the slightly different model, which includes 
evaporation, by Newman and Grima [34] and Grima |13j . 



B. Negative autochemotaxis 

First, we consider the microorganism's motion with no 
fluctuations present. In this case, the swimmer reaches 
a steady state at infinite time, which is determined by a 
constant swimming speed Vg{t — s- oo) = v^, where the in- 
dex s denotes the steady-state configuration. Under this 
condition, the drag force induced by the solvent, 7Vs, 
equals the driving force due to the chemical, — AVp. Fig- 
uratively, the microorganism surfs down its own chemore- 
pellent, which it excreted at times t' < t — to. Using 
Eq. ^ , the steady-state velocity Vs is determined by the 
self-consistent equation 



7^^^ = 



2D, 
47ri/2n 



dt 



2A 



,exp [-v^,t'/{AD,)] 



{ATTDct'Y/'^ 



r 1- 



2' 4D, 



(16) 



where T{a, z) denotes the incomplete gamma function. 
Eq. ( 16 1 has non-zero solutions w,, > for any value of 



A in one and two spatial dimensions, whereas there is 
a dynamical phase transition at a lower critical value of 
A* = — Stt'^/'^ZJc^^^o^^T in three dimensions, below which 
the microbe comes to a rest. The latter is obtained an- 



alytically by expanding Eq. (16 1 up to second order in 
Va- For the parameters used in the simulations reported 
above, = 10, 7 = 10, and to — 0.001, the transi- 
tion occurs at A* « —4450. The asymptotic solutions for 
small |A| in one and two dimensions or for small |A — A*| 
in three dimensions, respectively, are given by 



(2i^c7)" |A| , 
2(i^,/to)'/'exp 



7e 
2 



d= 3, 



wher 



(17) 



e 7_B 



2.7183 is Eulcr's constant. The steady-state 
velocity as a function of |A| is shown in Fig. |4] for all 
three spatial dimensions. For high coupling constants. 




FIG. 4: (color online). The steady-state velocity Vs{\) as a 
function of coupling constant |A| for one, two, and three spa- 
tial dimensions, determined theoretically (lines), compared to 
the square root of the slopes of the mean-square displacements 
at early times, dyj {{rb{t) — r(,(0))2) /dt, for two and three di- 
mensions (symbols) [cf. Fig. [2](b, c)]. Inset: The asymptotic 
solution of lis for small values of | A — A* | , where we set A* = 
for d= 1,2. 



A > 10**, the steady-state velocities Vg agree well with 
the square root of the slope of the mean-square displace- 
ments Vs ~ (3a/ ((rh(t) — rt,{0))'^)/dt, as obtained from 
the simulations with noise in two and three dimensions 
(see also Fig. ^ . 

If noise is put on, again, the microbe is eventually dis- 
turbed in its steady-state motion. In one dimension, the 
microbe needs to overcome a barrier in order to change 
the direction of motion from left to right or vice versa, 
which renders the steady state very stable. As we were 
not able to observe long-time diffusive motion in one di- 
mension for any coupling constant but always found bal- 
listic motion within the time window accessible in com- 
puter simulations, we do not attempt to give a theoret- 
ical estimate of the diffusivity. In two or three dimen- 
sions, the picture is very different: The microorganism is 
only constrained in its motion parallel to its current tra- 
jectory, whereas it is free to move perpendicular to the 
same. During such transverse fluctuations the direction 
of the microbe's motion changes with the gradient of the 
chemical fleld. 

In the following, we give a theoretical estimate of the 
microbe's long-time diffusivity int c? = 2,3 under the as- 
sumption of a fast relaxation of the chemical fleld, which 
is justifled for ^ D. In this case, the change of 
orientation is determined by the local, time-independent 
curvature kqW of the isodensity line (in 2D) or surface 
(in 3D) at the steady-state position rs{t); the isodensity 
planes are deflned by {r*(i)| p(r*(t),t) ^ p{rs{t),t)}. Lo- 
cally, they have the form of a parabola (in 2D) or of a 
paraboloid of revolution (in 3D) and they move with the 
microorganism at the microorganism's velocity; this can 
be seen in Fig. |5] where a typical trajectory and the ac- 
cording chemical fleld at time t is plotted for a microor- 
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FIG. 5: (color online), (a) Snapshot of the instantaneous 
density profile p{x,y) of the chemorepellent released by the 
microorganism moving in 2— dimensions, obtained from sim- 
ulation, at time instant t — 10.0 (in units of Ae~^). (b) The 
entire trajectory, shown as the red (thick) curve, of the mi- 
croorganism. The current position of the microorganism r(,(t) 
is indicated by the blue (black) dot in both the figures, and 
the direction of motion is indicated by arrows along the trajec- 
tory. The corresponding coupling strength being |A| — 10000. 



ganism in two dimensions. 

Fluctuations transverse to the direction of motion lead 
to a mean transverse displacement {r']^{t)) — 2{d~ l)Dt, 
which, in turn, leads to a change in orientation of the 
velocity director v;, by a mean-square angle 



{0\t))^4{vlit))^2id-l)D4t. 



(18) 



Exploiting that 6 is Gaussian distributed, the average, 
normalized, and on the initial orientation Vb(t = 0) pro- 
jected velocity vector Vb(t) is given by 



(Mt) ■ MO)) 



= (cos[6'(t)]) = exp 



{0'{t)) 



(19) 



Replacing time by arc length s — Vgt, the problem un- 
der study can therefore be mapped to the worm-like 
chain (WLC) model [H [THl [35] of a polymer of a to- 
tal arc length L = Vgt and a persistence length, which 
is defined implicitly by (cos[0(t)]) = exp[— w^i/Zp]. Mak- 
ing use of Eqs. (18) and (19 1 the latter therefore reads 
= ^s/[('^ ~ 1)I?Kq]. The mean square end-to-end 
distance of the WLC for long chains L ^ Ip is well 
known ^ to be given by 



{4it)) 



'2Llp — 



2vit 



{d~l)DKl 



(20) 



The long-time diffusivity of the microbe is then obtained 
as the time derivative of the mean square end-to-end dis- 
tance. 



A (A) 



1 d{vl{t)) 
2d dt 



d{d-l)DKl{\) 



(21) 



where we point out the A-dependence of the steady-state 
velocity and the steady-state curvature. Di{\) is plotted 
m Fig. [2] (d) and compared to the results of the simula- 
tion. 

It is ascertained from Fig. |2] (d) that, in the limit of 
large A > 10^, the theory overestimates the diffusivity by 
a factor of ^ 3, both in two and three dimensions, for the 
parameters D = IQ and to = 0.001, almost independent 
of A. This discrepancy is attributed to the strong as- 
sumption of a constant, non-fluctuating curvature ko(A) 
in Eqs. ( 18 ) and ( 21 1 . Relaxing this assumption to fluctu- 
ating curvature with first and second moments (k) = Kg 
and ((k— Ko)^) — Ak^, which is further uncorrelated with 
the transverse displacement, i.e., {K{t)r±{t')) = 0, the 
diffusion constant reduces to Di -D;/[l + (Ak/kq)^]. 
Therefore, although we do not attempt to give an esti- 
mate of Ak here, the long-time diffusivity is expected to 
be smaller than in the "zeroth order" theory, Eq. (21 1, 
in agreement with the simulation results. 



V. CONCLUSION 

In conclusion, we have explored the dynamics of au- 
tochemotaxis: a model microorganism is "smelling" its 
own secretion which is diffusing away. The microorgan- 
ism follows the gradient [30] of its secreted chemical. 
For the attractive case, the mean-square displacement 
of the microorganism reveals a transient dynamical ar- 
rest, most pronounced in low spatial dimensionality. In 
the opposite case of chemorepulsion, there is a transient 
ballistic behavior which crosses over to ultimate diffusion 
[37} [38] . A simple theoretical analysis for large coupling 
strengths A reveals a scaling law for the long-time diffu- 
sion in the case of a chemoattractant, by regarding small 
excursions of the microbe's position about the bottom 
of an effective, time-dependent trapping potential. In 
the case of chemorepulsion, the microbe's trajectory was 
mapped to the contour of a worm-like chain, which gives 
a semiquantitative agreement with our computer simu- 
lations. The ranges of the coupling strength reported 
here are also easily obtainable in real experimental sit- 
uations (A ~ 10~^ — lO"*), as estimated earlier in the 
text. We note that though the mechanism of temporal 
sensing [22l|32j|33] of chemoattractant or chemorepellent 
which we do not consider here is found in many small, 
fast moving bacteria like E. coli, the direct gradient sens- 
ing mechanism is also present in microorganisms like the 
amoeba Dictyostelium discoideum, the yeast cell Saccha- 
romyces carevisiae, lymphocytes, glial cells and myxobac- 
teria jMllSnilSI]- 
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Chemotaxis in the gliding bacterium M. xanthus was 
demonstrated to be in response to self-generated sig- 
naling chemicals [39] • However, the difFusivity of the 
chemoattractant is significantly lower compared with 
the bacterial motility in this peculiar case. A numer- 
ical investigation |:40^ of a simple model case in one- 
dimension compared the interplay between chemotaxis 
and chemokinesis mechanism with a concentration de- 
pendent switching rate, showing crossover from sup- 
pressed to enhanced diffusion in the parameter space. 
This mean-field approach, studied in the limit of vanish- 
ing chemical diffusivity and chemical degradation rate, 
was based on the simplifying assumption that the chemi- 
cal coupling only affects the frequency of direction rever- 
sal of the cell, keeping the cell speed unaltered. 

Models of active colloids using surface reactions as a 
potential mechanism for self-propulsion has been pro- 
posed (see e.g. [4Hl42|). In Ref. [4r, the model molecular 
machine is a spherical colloid which reacts with the sub- 
strate at a specific site on its surface and self-propel by 
asymmetric release of the reaction product. The time 
evolution of the product particles is similar to that of 
chemical molecules emitted by a microorganism. Though 
the exact way of imparting a biased motion to the col- 
loid by the emitted particles is very different from chemo- 
tactic coupling, such systems may also provide interest- 



ing situations to compare with our results. In fact in 
Ref. ^42j , active colloids interacting with a self-generated 
cloud of solute were shown to have distinct propulsive 
and anomalous super-diffusive regimes preceding a final 
long-time effective diffusion. This is worth comparing 
with the ballistic to a modified diffusion crossover in our 
study on negative autochemotaxis. 

Future work should focus on generalizing the model 
to a collection of model microorganisms where 
hydrodynamic flow effects can play an important role 
[23, 43, 44 and steric repulsions can lead to aggregation 
and clumps [J^ as known from active particles [24l |46] . 
It would further be interesting to study the case of a ge- 
ometric confinement of the particles |371 HHl SS] and the 
chemoattractant/repellent in order to see effects of a di- 
mensional crossover. Finally, improving the theoretical 
approaches towards full quantitative agreement will be 
an important task for future research. 
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